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Signatures of magnetar central engines in short GRB lightcurves 
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ABSTRACT 

A significant fraction of the Long Gamma-ray Bursts (LGRBs) in the Swift sample have a 
plateau phase showing evidence of ongoing energy injection. We suggest that many Short 
Gamma-ray Bursts (SGRBs) detected by the Swift satellite also show evidence of energy in- 
jection. Explaining this observation within the typical SGRB progenitor model is challenging 
as late time accretion, often used to explain plateaus in LGRBs, is likely to be absent from the 
SGRB population. Alternatively, it is predicted that the remnant of NS-NS mergers may not 
collapse immediately to a BH (or even collapse at all), forming instead an unstable millisec- 
ond pulsar (magnetar) which powers a plateau phase in the X-ray lightcurve. 

By fitting the magnetar model to all of the Swift SGRBs observed until May 2012, we find 
that about half can be clearly fitted with a magnetar plateau phase while the rest are consistent 
with forming a magnetar but the data are insufficient to prove a plateau phase. More data, both 
at early times and a larger sample, are required to confirm this. This model can be tested by 
detecting the gravitational wave emission from events using the next generation gravitational 
wave observatories. 
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1 INTRODUCTION 



Following the launch of the Swift satellite dOehrels etaDl2004l) . 
it has been possible to place tighter constraints on the na- 
ture of short gamma-ray bursts (SGRBs). The detection of their 
faint and rapidly fading X-ray afterglows has led to the iden- 
tification of optical afterglows and in ma ny cases candidate 
host galaxies (for example GRB 050509B, lOehrels etal1l2005l : 
lHiorthetal.ll2005h . These observations have provided significant 
support for the popular compact binary merger progenitor theo- 
ries, i.e. the coalescenc e of two neutron stars (NS) or a NS and a 
black hole (BH) (Lattimer & Schrammlll976l : lEichler et al.lll989l : 
iNaravan. Paczvnski. & PiranI Il992h . However, without the coin- 
cident observation of gravitational waves by observatories like 
LIGO (Laser Interferometry Gravitational-wave Observatory) we 
are missing the supporting "smoking gun" observation for this pro- 
genitor theory. 

Observed features in X-ray lightcurves suggest longevity of 
the central engine of GRBs, for example late time flares (e.g. 
[Curran et al. 2008; Margutti et al. 2010; Bemardini et al. 201 1) and 
plateaus (e.g. Nousek et al.l2006l : IZhang et al.l2006h . GRBs whose 
X-ray lightcurves have a steep decay and a plateau phase followed 
by a standard afterglow phase, have been identified as "canonical" 
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lightcurves (iNou sek et al.' '20061 : lO'Brien etaLll2006l : IZhang et al.l 
l2006l : lEvans et al.li200 9). The steep decay phase is associated with 
high latitude emission from the prompt emission f ollowed by a 
late emissio n plateau giving the plateau phase ( Tagliaferri et al.l 
l2005l : iGoad et al. 2006) . The fluence of this plateau can be com- 
parable to the fluenc e of the prompt emission (lO'Brien et al.ll200^ : 
lMarguttiet"aDl2012l) . and typically they occur from 10^ - 10^ s 
till 10^ - 10^ s after the trigger time. The plateau is thought to 
provide evidence of ong o ing central engin e activity (iNousek et al.l 
l2006l : IZhang et al.ll2006l) . lEvans et all (12009 ) studied 162 GRBs in 
the Swift sample identifying a "canonical" lightcurve in 42% of 
GRB X-ray Ughtcurves, including 2 (051221A and 060313) out of 
11 SGRBs analysed. 

Although studies of flares and plateaus are typically con- 
ducted for LGRBs, fainter versions are evident in many SGRB 
X-ray lightcurves s uggesting a long lived central engine (e.g. 
iMargutti et al.l201ll) . This is problematic for SGRB progenitor the- 
ories as accretion is expected to end within a few seconds (pow- 
ering relativistic jets; Rezzolla et al. 2011) and only a small frac- 
tion of the merger mass is available (0.01 - 0.1 M(;) although this 
is dep endant on the NS equation of state, iLee & Ramirez-Ruia 
l2007h . Additionally, it is thoug ht that the accretion disk ^ets de - 
stroyed after a few seconds (e.g. lMetzger. Piro. & Ou ataert 2008). 
There have been studies of fallback accretion, in which the NS is 
shredded and parts (^ 10% of the original disk mass) are flung 
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Table 1. The Swift SGRB sample and the results of broken powerlaw fits to the observed BAT-XRT data in the 0.3-10 keV band (as described in the text) and 
the X-ray spectral indicies for each regime (Vx). These are subdivided into those with 2 or more significant breaks in their lightcurves, those with 1 break and 
those with no significant breaks. Where values are left blank there was insufficient data available to constrain them. The last column shows the probability that 
this fit is a chance improvement on a simpler model. 



into highly eccentric orbits which accrete onto the c entral en- 
gine at late times giving flares in the X-ray lightcurve (iRosswod 
12007 '). Flares may also be caused by Toomre insta bilities within 
the accretion disk IPema. Armitage. & ZhandEo06h . although this 
does not explain plateau emission or late time flares as the SGRB 
accretion disks are expected to accrete within the first few sec- 



onds . ICannizzo. Troia. & Gehrelsl(EoT l|) have attempted to explain 
plateaus by introd ucing a band of material at a large d istance from 
the central engine. 'C annizzo. Troia. & Gehrelsl (I2OIII) suggest that 
the required reservoir of material could be provided via the accre- 
tion disk moving outwards (due to having a large amount of angu- 
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Table 2. Properties of the SGRB sample, including Tgo, and Fluence (15-150 keV). These observed quantities, including host galaxy associations, offsets 
and optical afterglow detections, are from published papers and GCNs (references listed below), host offsets are quoted with errors if published. When the 
redshift is not known, the average redshift 0.72 was used and this is shown using brackets. 



(i)Cummings et al. (2005b); Soderberg et al. (2006) (2)"Markwardt et al."( 2006a") ;"Roming et al."("2006) (3)'Markwardt et al." ('2006b') ;"Stratta et al."('2007|) 
(4)Parsons et al. (2007); Berger et al. (2009); Kocevski et al. (2010) (5) Sato et al. (2009); AntonelH et al. (2009); Xin et al. (2011) (6)Barthelmv et al] 

( 2009); Rowlinson et al. ( 2010a) ('^)Barthelmv et al. (2010a); Tanvir & Levan (2010) (8)Baumgartner et al. (2010) (9)Krimm et al. ( 2010^^ 

IChomock & Berger (201 1) (lO)Palmeret al. (2012) (ii)Cummings et al. (2012) (i2)Sato et al. (2005b); La Parola et al. (2006) (i3)Sato et al. (2006^ 
Cucchiara et al.l (l2006h ( i^)Barthelmv et al. (2007) (i5)Krimm et al. ( 2007); Perlev et al. ( 2008) (i6)Cummings et al. ( 2008a) (i'^)Cummings et al. ( 20^ 
Rowlinson et^ 1201 Obi) (i8) .Baumgartner et al. (^20081) : [Immler & Holland (.2008.):.Covino et al. t2008>) (iQ).Uk watta et al. ( 2009); de Pasau ale et al..aoToll 
iMcBreen et al.l ( l2010[) (^o ferimm et al.l ( l2009b[) (^ ^Markwardt et al.l (l2009h : Levan et alT(l2009h : lMales^ e t al. ( 2009) (^2) Sakamoto et aklhoi lbl) 



lar momentum) or ejecta thrown out during the merger in highly 
eccentric orbits that circularises forming an accretion disk. 

An alternative theory is that during some GRBs a millisecond 
pulsar (magnetar) may be formed with eno ugh rotati onal energy 
to prevent gravitational c ollapse (,Usov. l992i: [Duncan & ThompsonI 



1 1991 lDai&Ljl998al lbl: IZhang & Meszarosll200lh . The rotational 
energy is released as gravitational waves and electromagnetic radi- 
ation, causing the magnetar to spin down. If the magnetar is suf- 
ficiently massive it may reach a critical point at which differential 
rotation is no longer able to support it, resulting in collapse to a 
BH. Assuming constant radiative efficiency, the energy injection 
from the magnetar would produce a plateau in the X-ray light curve 
Jzhang & Meszaros 2001) and would be followed by a steep decay 
if the magnetar collapses to a BH. The progenitor of this system is 
typically tho ught to be a collaps ar a nd LGRB candidate s have been 
identified bv lTroia et all (l2007h and lLvons et~aD (l201Qh . However, 
it has also been proposed that such a magnetar could be formed 
by the merger of two ne utron stars (Dai & Lu 1998a; Dai et al. 
l2006l : lYu& Huand 120071) or via the accretion induced collapse 
(AIC) of a white dwarf (WD) dNomoto & Ko ndo" 1 99 ll: IUsovl 19921 : 
iLevan et alJl2006l : lMetzger. Ouataert. & Tho mpson 2008|). A can- 
didate event for this is GRB 090515 with an unusual X-ray plateau 
followed by a steep decay (Rowlinson et al. 2010a). The likeli- 
hood of producing this event is dependent on the equation of state 



of neutron stars. iMorrison. Baumgarte. & Shapirol (|2004|) studied 
the effect that the equation of state of a NS and rotation would 
have on the remnant of a compact merger, i.e. whether a NS or a 
BH is formed (see also Shibata & Taniguchi 2006). They showed 
that, even for the harder nuclear equations of state, the rotation of 
the NS could increase the maximum mass by ^ 50% and hence 
mergers could often result in a NS. Considering the parameters of 
6 known Galactic NS binaries and a range of equations of state, 
Morrison, Baumgarte, & Shapiro (2004) predict that the majority 
of mergers of the known binaries will form a NS. 



^The recent discovery of an 1.97 M© NS jDemorest et alj 

'2010^ provides further supporting evidence of the possibility that 
high mass magnetars can be formed from NS mergers (the max- 
imum mass of NSs is dependent on the very uncertain NS equa- 
tion of state, s o this is a conservat ive lower limit on the maximum 
mass of a NS). IOzel et alj (l2010bh show that, for a maximum non- 
rotating NS mass of Mmax =2.1 M©, the merger of two NSs with 
a total mass ^ 1.4Mmax will have a delayed collapse to a BH 
(i.e. a magnetar phase). They also predict a regime in which the 
merged remnant does not collapse to form a BH, in this case the 
total mass is ^ 1.2Mmax. If the maximum NS mass is 2.1 M©, 
then the merger of two NSs of masses up to 1.3 M© would result 
in a stable magnetar and the merger of two NSs with larger masses 
(up to 1.5-1.7 M©) would form an unstable magnetar. As the ma- 
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(0.72) 


0.12±0.03 


0.63±0.17 


1.4±0.2 


7 




N 


100628A(42) 


(0.72) 


0.36±0.009 


1.26±0.25 


0.3±0.1 


y 




N 


110112A(43) 


(0.72) 


0.5±0.1 


2.14±0.46 


0.3±0.1 


y 




Y 


111117A(44) 


(0.72) 


0.47±0.09 


0.65±0.22 


1.4±0.2 


y 


1.00±0.13 


N 



Table 2 - continued 



(24) |Barthelmv et al.l ( l2005bl):lGehrels et al.l (|2005|) (25)|Sato et al.| |2005j) (26)|Parsons et al.](l2005h:lLevan et al.l (l2008h, note this is a candidate extra-galac tic 
magnetar giant flare (2'^)Cummings et al. ( 2005a); Barbier et al. (2005); Klose, Laux, & Stecklum ( 2005) (28) Sato et al. ( 2006a) (29)"Parsons etal."(200^; 
Ziaeepour et al. (2006) (3°) Sakamoto et al. ( 2007a) (3i)Tueller et al. (2007); Antonelli et al. (2007); HoUand, de Pasquale, & Markwardt (2007) 
(32)|Sato et al. ( 2007a); Berger & Kaplan ( 2007) (33) Sakamoto et al. ( 2007b); T hoene et alt ( 2007) (34)Fenimore et al. ( 2007) (35)Krimm et al. ( 200^ 

(36) Barthelmy et al. ( 2008) (3'^) Barthelmy et al. ( 20Q8b) (38 ) Automated BAT analysis products (39 ) Krimm et al.. i2009a); Cenko et al. ( 2009) 

(40) |Markwardt et_£.l (l2010l) jLevan et al.l(l2010l):lFone 61^111(120101) (^i) Sakamoto et al. ( 2010); Miller et al.l(l201Ql):lLevan et al.l (l20ld):lBerger et al l(l2010l) 



(^2) |Bartlielmv et al.l ( l2010bh;IStarling. Beardmore, & Immleil (1201 Oh;lBerger (2010.) (^3) Barthelmy, Sakamoto & StamatikosI |201lh ; lLevan. Tanvir & Baked 
1201 ih (^^) |Sakamoto et al.l i2011al) : ICenko & Cucchiaral ( l201lh ; lBerger. Pong. & Sakamoto! i201lh 



jority of observed NSs have masses ^1.4 M©, it seems reasonable 
to predict that many NS mergers could result in a magnetar. The 
stability of the final magnetar is dependent on the maximum possi- 
ble mass of a NS which is still uncertain. Its lifetime depends both 
on the rate that additional mass (if any) is accreted after formation, 
as well as the rate at which angular momentu m is extracted by e.g. 
gravi t ational waves or magnetic torque s (e.g. IShibata & Taniguchil 
l2006l :IOechslin. Janka. & Marek 2007). 

In this paper, we consider all Swift detected SGRBs, Tgo ^ 2 
s, observed until May 2012 with an X-ray afterglow or which were 
promptly slewed to and observed by the X-ray Telescope (XRT; 
iBurrows et al.ll2005l) . This allows the inclusion of SGRBs without 
an X-ray afterglow but which do have a constraining upper limit. 
For all the SGRBs, we analysed the BAT (Burst Alert Telescope; 
iBarthelmv et al .12005 ah data by creating lightcurves using a variety 
of binning in signal-to-noise ratios and time looking for evidence 
of extended emission at the 3cr level where we consistently saw ex- 
tended emission over more than 30 s (the SGRBs with identified 
extended emission are 050724, 050911, 051227, 060614, 061006, 
061210, 070714B, 071227, 080123, 080503, 090531B, 090715A, 
090916 and 111121 A). This procedure recovered all of the ex- 
tended emission bursts identified by Norris, Gehrels, & Scargle 
(|201Q|) . Hence, our selection criteria excludes SGRBs with ex- 
tended emission, which may share a common progenitor to SGRBs 
but this remains uncertain. This sample is used to identify those 
with a plateau phase in their lightcurves suggesting ongoing cen- 
tral engine activity. These results are discussed in section 2. A sub- 
sample with sufficient data are then studied for the signature of a 



magnetar (with or without collapse to a BH) which may signify 
the coalescence of two NSs. If found, this would provide additional 
support to this popular progenitor theory although forming a mag- 
netar via the AIC of a WD is not ruled out. The magnetar model is 
considered in section 3; with a description of the model and sample 
used and analysis of the available data. A discussion of the im- 
plications, e.g. for gravitational waves, is given in section 4 and 
our conclusions are given in section 5. Throughout this work, we 
adopt a cosmology with Hq — 71 kms~^ Mpc~^, Qm = 0.27, 
^^A = 0.73. Errors are quoted at 90% confidence for X-ray data 
and at Icr for fits to the magnetar model. 



2 PLATEAU PHASES IN SGRB LIGHTCURVES 

Out of our sample of 43 SGRBs, shown in Table [T] only 6 did not 
have a detected X-ray afterglow (GRBs 050906, 051105, 070209, 
070810B 071112B and 081101). Hence, -86% of Swift SGRBs 
with a prompt slew have detected X-ray afterglows. The observed 
properties of the SGRB sample are given in Table |2] 

The 0.3 - 10 keV observed flux X-ray lightcurves were ob- 
tained from the automated analysis page for eac h individual SGRB 
from the UK Swift Science Data Centre website (lEvans et al .120071 . 
l2009l) . The Burst Alert Telescope (BAT) lightcurves were created 
using standard pipelines in the He AS OFT package with 3cr signifi- 
cance bins. The 15-150 keV BAT spectra were fitted in XSpec for 
each SGRB and then extrapolated to obtain the flux at 0.3 - 10 keV. 
Using this extrapolated flux and the net count-rate in the BAT spec- 
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Figure 1. These are the BAT-XRT lightcurves (0.3 - 10 keV, observed flux) sorted into 3 groups, a) These GRBs have 2 or more breaks in their lightcurve. b) 
GRBs with 1 break in their lightcurve. c) and d) GRBs with no significant breaks in their lightcurve. 



trum, each count-rate data point in the BAT lightcurve was scaled 
to an 0.3 - 10 keV flux using a simple power law spectral model. 
These were combined with the XRT lightcurves to make the BAT- 
XRT lightcurves used in this analysis. For later comparison with 
the magnetar model, these lightcurves were converted into unab- 
sorbed flux lightcurves and then into re stframe 1-10000 keV lumi - 
nosity lightcurves using a k-correction dBloom. Frail. & Sa rj200lh 
giving an approximation to a bolometric lightcurve. The range of 
k-corrections obtained are typically consistent with the range 0.4- 
7 as obtained by Bloom, Frai L & S aril (1200 It) . However, there are a 
small number of large k-corrections (particularly for GRBs 070809, 
080905 A and 1012 19A), showing the spectrum may be poorly un- 
derstood at the frequency range we extrapolate to. If no redshift 
is known the mean SGRB redshift is used, z ^ 0.72 (excluding 
the redshift for GRB 061201 as t he host galaxy association re- 
mains uncertain: IStratta et al Jl2007l Tunnicliffe et al. in prep), and 
the implications of choosing this average redshift are discussed in 
Section 3.3. All the SGRB observed BAT-XRT lightcurves were 
fitted with multiple power laws from the final decay phase in the 
BAT prompt emission throughout the total X-ray afterglow using 



These fits were then used to identify those with a "canoni- 
cal" like lightcurve. An XRT spectrum was created for each region 
of the lightcurve using the aut omatic data products on the UK Swift 
Science Data Centre website (lEvans et al. 2007, 2009). The SGRB 
lightcurves are shown in Figure [T] 

We assume that Fj^ oc iy~^t~^ where /3 = F — 1 is the spec- 
tral index, F is the photon index (F^ is the photon index measured 
using BAT and F^c is the photon index measured using XRT) and 
a is the temporal index. We define the steep decay phase following 
the prompt emission to have a power law decay of ai, after which 
the decay can break to a decay of a2 and a further break to as- We 
always define a2 to be the shallowest decay phase, as this allows 
the direct comparison of all plateau phases in the subsequent anal- 
ysis. All SGRBs with 1 or more breaks in their lightcurves have 
a plateau phase. In three cases there are more than two breaks in 
the lightcurve. GRB 070724 has a third break at T3 = 152^5^ 
(a4 = l.lSlo.Se and F4 = ^-^^to.tll GRB 090515 has a third 
break at T4 = 24lt^o (0^4 = lOto.g?) and GRB 101219A has a 
break at T3 = 241^^^ (0^3 = l-88lo.25 and F4 = 1.63t° j^). 

^ https://heasarc.gsfc.nasa.gov/docs/software/ftools/others/qdp/qdp.html 
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Figure 2. Histograms showing the break times for the SGRB lightcurves 
with a plateau phase. Ti is the break from the steep decay phase to the 
plateau phase while T2 marks the end of the plateau. The blue filled his- 
tograms correspond to the SGRB sample u sed in this Paper an d overplotted 
in red are the LORE values determined bv lEvans et al .1 ( 1200^ . 



Figure 3. Histograms showing the temporal indices of the SGRB 
lightcurves with a plateau phase, ai is the initial steep decay phase from 
the last decay in the prompt emission. 0:2 are the plateau and shallow decay 
phase slopes. 0:3 is the final afterglow decay slope. The filled histograms 
correspond to the SGRB sample us ed in this pa per and overplotted are the 
LGRB values determined bv iEvans et all ( l2009h . 



The 6 GRBs which were undetected by XRT are fitted with lower 
limits for ai using the shallowest decay allowed by the BAT data 
and the XRT upper limit. In Table [T] we provide the lightcurve fits 
for all the SGRBs in the sample. An F-test was conducted using 
the ^nd degrees of freedom for each fit to determine the proba- 
bility that the fit is a chance improvement on a simpler model (i.e. 
an F-test between the model provided in Tableland a model with 
1 less break in the lightcurve). We utilise the method described in 
lEvans et aP (l2009h to determine the best fit, i.e. the model that has 
the most breaks and the probability of being a chance improvement 
on a simpler model as ^0.3%. 

There are several caveats which need to be considered with the 
results in this Section and for the magnetar fits in Section 5.3.2. As 
SGRB afterglows are often faint and fade rapidly, these lightcurves 
and spectra can be poorly sampled giving large errors on the val- 
ues in Table [T] This could also cause breaks in the lightcurve to 
be missed due to large bin sizes (bins typically contai n 20 photons 
in PC mode data so bins could have long durations; lEvans et aP 
(2007). Additionally, the Swift satellite slews to observe GRBs af- 
ter detection, leading to a characteristic gap between the BAT data 
and the XRT data, and XRT can only observe for short windows. 



due to Earth occultation during orbits, giving further gaps in the 
lightcurves which could also hide features in the lightcurves. 

Using the broken powerlaw fit method, we find that 22 
SGRBs (^50%) are consistent with having a plateau phase in their 
lightcurves, although the plateau phase is not always directly ob- 
served due to the gap in the lightcurve prior to the XRT observa- 
tions. It is hard to rule out plateau phases in other cases (since the 
plateau phase could be missed by the sampling or lost due to the 
faintness of the afterglow). Those which were undetected by XRT 
do not require extreme decay slopes relative to the rest of the sam- 
ple of SGRBs. The break times of the SGRBs with plateaus are typ- 
ically occuring orders of magnitude earlier than for the canonical 
LGRBs (as shown in Figure [2]). This may be ca used by our use of 
BAT and XRT data whereas lEvans et all (l2009h use only the XRT 
data and is only able to find plateaus at times after XRT has started 
observing. However, it is very rare for XRT to not observe the steep 
decay phase for LGR Bs so the inclusion of BAT data does not affec t 
the plateau fits (e.g. (o^rien et all 1200^ : IWillingale et al.ll2007h . 
Additionallv JEvans et al.l (l2009h discussed whether their type b and 
type c LGRBs can be cannonical (i.e. those which are steep then 
shallow or shallow then steep). They conclude that they are not 
based on the plateau decay rates and break times for type b and the 
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Figure 4. (a) The BAT fluence (15-150 keV) plotted against the XRT un- 
absorbed flux at 100 s (0.3 - 10 keV). Blue stars have 2 or more significant 
breaks in their lightcurves, green circles have 1 break and red triangles have 
no significant breaks in their lightcurves. (b) The BAT peak photon flux (15- 
150 keV) against the XRT unabsorbed flux at 100 s in the observer frame 
(0.3-10 keV). Symbols are as in (a). 



relative BAT versus XRT fluxes for type c. Whereas for SGRBs, the 
BAT observations need including in order to identify the steep de- 
cay phase. Histograms showing the various SGRB decay slopes for 
those with plateau phases are sho wn in Figure [3] with t he values for 
canonical LGRBs, deteraiined bv lEvans et alJ (l2009l) . The values 
for ai and 0^2 are consistent with the LGRB sample, but the final 
decay phase (0^3) is typically steeper than fo r the LGRB counter- 
parts (consistent with the results obtained by iMargutti et alj|2012l . 
where the overall decay of SGRBs is typically found to be steeper 
than LGRBs, however their sample of SGRBs used for this is dom- 
inated by SGRBs with extended emission). Using a Kolmogorov- 
Smirnov test between the values for LGRBs and SGRBs, ai is 
consistent with being drawn from the same distribution (p-value = 
0.07), although the values for 0^2 are unlikely to be from the same 
distribution (p-value = 0.003) and 0^3 are highly likely to be drawn 
from completely different distributions (p-value = 0.00007). In the 
following analysis we consider SGRBs with 2 or more breaks in 




Prompt Fluence {15 ■ 



10^ 

1 50 keV) (erg cm-^) 



Figure 5. The prompt BAT 15 - 150 keV fluence in comparison to the 
shallow decay phase unabsorbed X-ray fluence extrapolated to the 15 - 150 
keV energy band. Symbols are as defined in Figure |4] and the black line 
shows where the shallow decay phase fluence is equal to the prompt fluence. 



their lightcurves (blue stars, the GRBs in Figure [1^), 1 break (green 
circles, the GRBs in Figure [TJ)) and those with no breaks in their 
lightcurves (red triangles, the GRBs in Figure [T]: and d). 

The BAT fluence (15-150 keV) of these GRBs is plotted 
against their 0.3 - 10 keV flux at 100 s in Figuregt- Those GRBs 
with a plateau tend to be clustered at somewhat higher fluence s 
and their X-ray fluxes are significantly higher at 100 s (^ 10~^^ 
- 10 ~^ erg cm~^ s~^). The GRBs which do not have a plateau 
phase in their lightcurves tend to have faint X-ray afterglows at 
100 s (^ 2 X 10~^^ erg cm~^ s~^) and relatively low fluences 
erg cm~^). Figure |4]d shows there is a wide variation 



(^ 2 X 10" 

in X RT flux at 100 s for SGRBs with similar prompt fluxes. 

lO'Brien et aP (|2006) and Willingale et al. (20(yi) found that 
the prompt fluence is comparable to the plateau fluence for LGRBs. 
In order to compare this result to our sample, we took the average 
flux for the plateau phase and multiplied it by the time at which 
the decay broke to a more typical afterglow (assuming this com- 
ponent started at the initial trigger time) giving the 0.3 - 10 keV 
fluence. This fluence was then converted to a 15 - 150 keV flu- 
ence using the spectral index and fitted absorption. Figure |5] shows 
the prompt and plateau fluence are generally comparable, which is 
consistent with the result obtained for LGRBs. There are four sig- 
nificant outliers (GRBs 061201, 070724A, 080905A and 090515), 
lying significantly above the one-to-one line, whose plateaus are 
significantly more energetic than their prompt emission. 

Figure [6] shows the spectral indicies plotted against the tem- 
poral indicies for the lightcurves with a plateau phase (values all 
tabulated in Table [B. Also plotted are the closure relations for the 
slow cooling reg ime (grey band) and the fast cooling regime (solid 
black lines) ( Zhang & Meszaros 2004). These show the same be- 
haviour identified by Evans et al. (2009) for the canonical sam- 
ple of LGRBs. In particular Figure [SJ) shows evidenc e of energy 
injection during the plateau phase as described by lEvans et aP 
(2009). These figures can be compared to updated values for the 
whole GRB sample using t he UK Swift Sci ence Data Centre 
(ww w.swift.ac.uk/x r t_live_ cat; lEvans et ID l2009h . 

IPainotti et al.l (l201Qh identified a correlation between the 
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Figure 6. The spectral index ^ versus the temporal index ol for the three 
regimes of the lightcurves with a plateau phase: (a) steep decay phase, (b) 
plateau phase and (c) standard afterglow phase. Where there is no XRT 
spectrum available for the steep decay phase, the BAT spectrum is used. All 
symbols are as defined in Figure [4] the solid lines and grey regions show 
the closure relations as defined bv lZhang & MeszarosI (l2004l) . and the black 
dashed line shows where o; = /3 + 2. 




1000 10- 
Restframe Plateau Duration (s) 

Figure 7. (a) The plateau phase unabsorbed flux versus the duration of this 
phase. Symbols are as defined in Figure|4] (b) The plateau phase luminosity, 
using published redshifts (filled symbols) or the average redshift (open sym- 
bols), v ersus the restframe du ration of this phase. The light grey data points 
are the Dainotti et al.l (|2pi0) sample of LGRBs. The black line shows the 
correlation between the luminosity and duration for the SGRB and LGRB 
samples which is consistent with the relationship found bv iDainotti et all 

( Eqiq.) . 



plateau phase luminosity and duration for LGRBs with a canonical 
lightcurve. Using redshifts where avai lable or the average SGRB 
redshift {z 0.72) and a k-correction dBloom. Frail. & SaribOQlh 
we calculated the luminosity and restframe durations for the SGRB 
sample (XRT fluxes used are the observed values which have not 
been corrected for absorption). These results are plotted in Figure 
|7]and the luminosity - duration correlation is identified. The fitted 
correlation for the SGRB and LGRB sample, 6 = -1.29 ± 0.12, 
log(a) = 48.74 zb 0.44, intrinsic scatter ay = 9 x 10"^^ =b 0.01 
(where Lx — a^piateau ^i^d the uncertainties on each datapoint 
and an intrin sic scatter are acco unted for in the fit using the method 
described in lD'Agostinill2005l) . is consistent with that for obtained 
the LGRB sample (-1.06 =b 0.28, 51.0 6 =b 1.02, although they 
did not account for the intrinsic scatter; D ainotti et al.lboioh . The 
SGRB plateau phases are typically more luminous and the plateau 
is shorter in duration than the LGRB counterparts. This may be a 
selection affect due to the inclusion of BAT observations in this 
analysis, and hence finding earlier plateaus. However, when BAT 
data are included in L GRB analysis the p l ateau properties do no t 
significantly change (lO'Brien et aLll2006l : Iwillingale et aDl2007h 
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additionally, there is a shortage of long duration plateaus observed 
in the SGRB sample. Cannizzo, Troia, & Gehrels (2011) argue that 
the relationship identified bv lDainotti et alj (|201Q|) is dominated by 
selection affects at z> 1.5 as there is an observational bias against 
faint plateaus due to the limiting XRT flux. However, SGRBs are 
typically at lower redshift (the SGRBs with an observed redshift in 
our sample have an average redshift of z^O.72) so our sample lies 
well within the region which is not dominated by selection affects. 

The plateau phases of GRB lightcurves are typically ex- 
plained as ongoing central engine activity, for example on go- 
ing accretion onto the central BH. However, ongoing accre- 
tion is problematic for NS-NS and NS-BH merger theories as 
there is insufficient surrounding material to maintain this accre- 
tion jLee & Ramirez-Rui^ l2007h . Fallback accretion from ma- 
terial on h ighly eccentric orbits has been postulated to re- 
solve this (IRosswo3 2007: iKumar. Naravan. & JohnsonI l2008l : 



ICannizzo. Troia. & Gehrelsl 2011 ). however it is unclear how to 
produce the required reservoir of material at a fixed radius. In the 
remainder of this paper, we suggest that the plateau phases could 
be powered by a magnetar formed via the merger of two NSs. 



3 MAGNETAR MODEL 

The magnetar model predicts a plateau phase in the X-ray 
lightcurve which is powered by the spin down of a newly formed 
magnetar. This section fits the model directly to the restframe 
SGRB lightcurves. The magnetar component is expected to be an 
extra component to the typical lightcurve. Therefore, we assume 
there is a single power law decay, ai, underlying the magnetar 
component. This value has been set to a i = F^ + l, where Tj is the 
photon index of the prompt emissio n, assuming that the decay s lope 
is governed by the curvature effect (iKumar & Panaitescul2000l) . i.e. 
that the surrounding medium is very low density as might be ex- 
pected for neutron star mergers. We note that the simple curvature 
effect assumed here does not account for any spectral evolution (for 
example as the peak energy moves through the observation band) 
however this is to keep the number of free parameters fitted in the 
model low. The normalisation of the power law decay fit is con- 
strained using the last decay from the prompt emission. In a small 
number of cases, the decay slope is significantly different from pre- 
diction and we allow ai to vary. It is important to note that the 
underlying lightcurve could be similar to other GRBs with a more 
complex afterglow light curve, but this work assumes that these are 
naked bursts (i.e. no surrounding ISM for neutron star mergers) and 
only the curvature effect is important. 

Also, we expect there to be flares also overlying the powerlaw 
decay and magnetar component (e.g. Margutti et al. 2011). Due to 
the limited statistics in SGRB lightcurves, we do not attempt to ex- 
clude possible flares from the lightcurve fits (except GRB 060313, 
which has multiple flares early in the lightcurve) and the underlying 
flares will slightly affect the fit parameters. 



3.1 Theory 

The model used here is as described in IZhang & MeszarosI fcOOll) 
and was s uggested to exp lain GRB 05 1221 A with a long lived 



and was s uggested to exp lain LrKB UMzzlA with a long lived 
magnetar dPan & X^?"2006'), for several L GRBs (iTroia et al.l l 20071 : 
iLvons et al.ll201Ql: Bemardin i et al.|[2012h and for the short GRB 
090515 jRowlinson et al. 2010a). This model is consistent with the 
late time residual spi n down phase driving a relativistic magnetar 
wind as described in iMetzger et all (120111) . We use the equations 



below with an underlying powerlaw component. Previously, the 
plateau duration and luminosity were calculated and then input into 
the equations. In this work, the equations are fit directly to the rest- 
frame light curves, taking into account the shape of the lightcurve 
(this is a comparable meth od to that used by IPalPOsso et al.l 
I2OIII : iBemardini et al.l 1201 2 who fitted a stable magnetar to the 
lightcurves of 4 LGRBs). We can then use the values of the mag- 
netic field and spin period obtained to derive the luminosity and 
plateau duration. 



Tem,3 — 2.05 (I45B ^^P? 



,15^0,-3^^6 
2 0-4 



Bp,15 — 4.2025/45/^6 ^/^0,49^em,J 
Po,-3 = 2. 05/45/^0, 49 ^em,^ 



(1) 

(2) 
(3) 
(4) 



where Tem,3 is the plateau duration in 10^ s, Lo,49 is the 
plateau luminosity in 10^^ erg s~^, I 45 is the moment of inertia in 
units of 10^^ g cm^, Bp, 15 is the magnetic field strength at the poles 
in units of 10^^ G, /^e is the radius of the neutron star in 10^ cm and 
Po,-3 is the initial period of the compact object in milliseconds. 
These equations apply to the electromagnetic dominated spin down 
regime, as the gravitational wave dominated regime would be ex- 
tremely rapid and produce a negligble electromagnetic signal. We 
have assumed that the emission is 100% efficient and isotropic as 
the beaming angle and emission mechanism remains very uncertain 
(see however section 3.4.4). The equations of vacuum dipole spin- 
down given above neglect the enhanced angular momentum losses 
due to neutrino-driven ma ss loss, which are ini portant at early times 
after the magnetar forms ("Metzger etal."201l'). Nevertheless, these 
expressions reasonably approximate the spin-down of very highly 
magnetized neutron stars of most relevance in this paper. Isotropic 
emission is also a reasonable assumption for relatively powerful 
magnetar winds, since (unlike following the collapse of a massive 
star) the magnetar outflow cannot be confined efficiently by the rel- 
atively small quantity of surroundi ng material expected following 
a NS merger or AIC ( Bucciantini et al.ll201lh . 

We use equation |5] to obtain the mass dependence of 
the model, where Mi 4 = IAMq, and equation [6] (from 
IZhang & Meszarosll200lh to determine the time dependence of the 
magnetar emission. 



/45 ^ MiaR{ 



Lem,49{T) — Lo,49 I 1 + 



10-3Te, 



(5) 
(6) 



If there is a steep decay phase after the plateau, it is assumed 
the magnetar has collapsed to a BH at the start of the steep decay 
(giving the Collapse Time parameter). The decay after collapse to 
a BH assumes the same powerlaw decay from the curvature effect, 
but starting at to = tcoiiapse- 

This model was then written into a QDP cod file (COmponent 
Definition file, used to generate new models within QDP which 
can then be fitted to data sets). In this analysis, the mass (Mi. 4) 
and radius (Rq) of the neutron star are constrained to be equal 
to 1 to reduce the number of free parameters in our model. These 
canonical values are consistent with the values determined by ob- 
servations of thr ee typical neutron stars, name ly M ^ 2M0 and 
7 ^ ^ 11 km dOzel. Bavm. & Guvejl2010ah . As the model con- 
siders an extreme neutron star, we note that the mass and radius 
may differ from these results. However, this only has a relatively 
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Figure 8. SGRB BAT-XRT restframe lightcurves fit with the magnetar model. The light grey data points have been excluded from the fit. The dashed line 
shows the power-law component and the dotted line shows the magnetar componenet. 
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Figure 8 - continued 
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Table 3. The SORB magnetar sam ple used with their magnetar fits. E^so, 1-10000 keV, is calculated using the fluences and redshifts in Table [J] a simple 
power law model and a k-correction iBloom. Frail. & SarilfcOO l). The values for ol are input into the model unless they are bracketed - in this case the values 
are fit within the model. If there is a steep decay phase, we assume the magnetar collapses to form a BH and the model determines the collapse time. The 
values for P_3 and B15 are fitted from the model assuming isotropic emission. Using the values of P_3 and B15 obtained from the model, we derive the 
plateau luminosity and duration using equations |2] and [T] The derived plateau duration is from the initial formation of the magnetar (i.e. the time of the GRB) 
and the point at which the X-ray emission from the magnetar starts to turn over from the plateau phase to a powerlaw decay phase. 



small aff ect on the mag netic fields and spin periods calculated (as 
shown in lRowlinson et al..2010a) and so it is a reasonable approxi- 
mation as we are just demontrating the plausibility of the magnetar 
model fitting the SGRB lightcurves. When this model is fit to the 
restframe lightcurves it produces Bp, 15, Po,-3, ol\ and the collapse 
time where appropriate. 



3.2 The Sample GRBs for magnetar fits 

The selected GRBs are those SGRBs in our sample with sufficient 
data to produce multiple data points in the X-ray lightcurve, giving 
a sample of 28 SGRBs. GRBs which have insufficient data to fit the 
magnetar model are not excluded from being magnetar candidates 
as it's possible to fit a range of realistic magnetar parameters with 
the minimal data points and unknown redshift. 68% of SGRBs in 
our samp le have been investigated for evidence of extended emis- 
sion by Norri s. Gehrels. & Scargle (2010) but, of these, none show 
evidence of extended emission. The remaining SGRBs in our sam- 
ple have no evidence of extended emission in their lightcurves (us- 



ing a variety of binning in signal-to-noise ratios and time looking 
for evidence of extended emission at the 3 a level). 

The magnetar sample are listed in Table|3] The restframe BAT- 
XRT lightcurves were fitted using the magnetar model, as shown 
in Figure [8] The lightcurves are fit over plateau region and the 
power law decay, including the last decay in the prompt emission 
and the X-ray observations. This removes the effect of the poorly 
understood flaring prompt emission not modeled by this method. 
We also provide the derived plateau luminosity and plateau dura- 
tion calculated using the magnetic field strengths, the spin periods 
and equations [1] and El The magnetar candidates fit the model well 
and the possible candidates are GRBs which may fit the magne- 
tar model if various assumptions are made. There are two potential 
outcomes: a stable long lived magnetar which does not collapse to 
form a BH and an unstable magnetar which collapses forming a BH 
after a short timescale (which have a collapse time in Table |3]). The 
following Sections compare the properties of the stable magnetars 
(blue stars in the figures) and the unstable magnetars which col- 
lapse to form a BH (green circles). We note that the fitted plateaus 
match the observations well but, due to insufficient data points par- 
ticularly prior to XRT observations, the plateaus are not always re- 
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quired by the observed data which can be fitted by simple broken 
power-law models. In some cases, the best fitting magnetar model 
gives a plateau phase ending prior to the start of the XRT observa- 
tions (e.g. 060801). In this situation, the fit is being constrained 
by the curving of the magnetar energy injection from a plateau 
phase to a powerlaw decline giving a characteristic curvature in the 
lightcurve (described by Equation 6). Therefore, the fitted model 
does not rely upon data during the plateau phase but instead uses 
the whole shape of the lightcurve. This leads to the model predition 
that those GRBs have a magnetar plateau phase which has not been 
directly observed, this can be used to test the model if we are able to 
observe SGRBs much sooner after the prompt emission with future 
X-ray telescopes. 

When fitting GRB 060 313, which may sho w evidence of late 
time central engine activity jRoming et al.ll200^ . it was noted that 
the model fits part of the lightcurve extremely well. In this case, 
we ignored the observations between 50 - 200 s (the initial X-ray 
data) in the fit as this duration appears to be dominated by flares. If 
these data are included in the fit, then the model does not fit the data 
well. The model fits well to GRB 090515 predicting values similar 
to those given in lRowlinson et al.l (l2010ah . 

In some cases, the model used here under predicts the flux at 
late times (for example GRBs 091109B, 100702A and 120305 A). 
This shows that our simple power law component, given by a sim- 
ple curvature effect model, is not sufficient and we should include 
spectral evolution or there may also be an additional afterglow com- 
ponent which has been neglected in this model. 



3.3 Analysis 

In Figure |9la) we show the spin periods and magnetic fields deter- 
mined for our sample of GRBs assuming isot ropic emiss i on. W e 
also plot the LGRB c andi dates identified bv ^Lvons et 
iDairOsso et'aD ( l201lh and lBemardini et al.l(l2012h . the SGRB can- 
didates tend to have higher magnetic field strengths and spin pe- 
riods. In Figure |9lb), we confirm the change in magnetic field 
strength and spin period caused by uncertainties in redshift ex- 
pected from previous analysis of GRB 090515 ( Rowlinson et al.l 
l2010ah . 18 of the SGRBs fitted by the magnetar model lie within 
the expected region of the magnetic field strength and spin peri- 
ods, these are the magnetar candidates listed in Table |2] 10 GRBs 
are outside the expected region (the possible candidates in Table 
12]). These GRBs may be in the expec ted (unshaded) region if they 
were at a higher redshift as shown in lRowlinson et al.l ( 2010ah and 
Figure |9lb). Additionally, this region is defined using angular mo- 
mentum conservation during the AIC of a WD (Usov 1992) and is 
not a physically forbidden region. Therefore, the candidates with 
spin periods >10 ms may remain good candidate magnetars. GRB 
051210 is included in the possible candidates list as it is spinning 
faster than is allowed in the models, but it is worth noting that if 
the NS formed had a mass of 2.IM0 then it would reside within 
the allowed region, as more massive NSs are able to spin at a faster 
rate. It is also worth noting that if GRB 051210 occurred at a lower 
redshift, as shown in Figure |9jb), or if the emission is significantly 
beamed then the spin period and magnetic field strengths would 
be higher and GRB 051210 would not be near to the spin break 
up period. The unstable magnetar candidates tend to have higher 
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Figure 9. (a) A graph showing the magnetic field and spin perio d of the magnetar fit s produced. The soHd (dashed) red fine and dark shaded area represent the 
spin break up perio d for a coUapsar (bin ary merger) progenitor |Lattimer & Prakashii2004 ) and the unshaded region sh ows the exp ected region for an unstable 
pulsar, as defined in lLvons e^^ |201 Q|) andlRowlinson et alj ( 1201 Oah . The initial rotation period needs to be ^10 ms ( IUsovlll992h and the lower limit for the 
magnetic field is G lThompsonl2007h . B lue stars = stable magnetar and Green circles = unstable magnetar wh ich collapses to form a BH. The black '+' 

symbols are the LGRB candidates identified bv lLvons et all fcOlOh i lPairOsso et"S1 l201lh : lBemardini et alJ ( l2012h . Filled symbols have observed redshifts 
whereas open symbols use the average SGRB redshift. (b) This graph is as (a) but focusing on the fits for two GRBs at different redshifts. The number below 
each data point is the corresponding redshift. GRB 060801 in blue is an unstable magnetar which collap ses to form a BH where as GRB 080702 A forms a 
stable magnetar. As expected, the paths of these lines are consistent with the predictions for GRB 090515 ( iRowlinson et al.ll2010ah . 



magnetic field strengths for their spin periods than the stable mag- 
netar candidates. The only exceptions are GRBs 1001 17A, which 
has been fitted with a stable magnetar model but would also be 
consistent with forming an unstable magnetar, and GRB 090426. 

3. 3. 1 Prompt and X-ray Properties 

In Figures [iota) and (b), the 0.3 - 10 keV flux at 1000 s or 10000s 
are compared to the flux at 100 s. The stable magnetar candidates 
tend to have a higher flux at 1000 s than the GRBs which are mod- 
elled as collapsing to a BH. This graph can be explained if we as- 
sume all SGRBs are occuring in a low density environment, re- 
sulting in little afterglow, and the only observed emission results 
from the curvature effect. The magnetar candidates which collapse 
to form a BH fade rapidly, whereas the stable magnetars are giving 
prolonged energy injection giving the higher late time X-ray fluxes. 
The stable magnetar candidate outlier in Figures [Tola) and (b) is 
GRB 1001 17A and it has already been noted that this GRB would 
also be fitted well by an unstable magnetar model. This analysis 
suggests that mergers collapsing straight to BHs have significantly 
fainter X-ray afterglows, which fade rapidly, and hence there may 
be a selection bias against these objects in our analysis (as we re- 
quired sufficient data points to fit the model). In Figures [TOlc) and 
(d) we plot the flux at 100 s and 1000 s versus the prompt 15-150 
keV fluence observed. At 100 s the unstable magnetar candidates 
clearly have a higher flux than comparable stable magnetar candi- 
dates (again GRB 100117A is the outlier) although this separation 
of the two populations has vanished by 1000 s. 

For each GRB in the sample, a 0.3 - 10 keV XRT spectrum 
(using th e automatic data produ cts on the UK Swift Data Centre 
website; lEvans et aLll2"007l . 120091) for the model derived rest frame 



plateau duration (converted to observed frame durations) was ex- 
tracted to compare the spectral properties in the proposed magnetar 
emission phase. This was not possible for some of the sample as 
XRT observations started after the plateau phase had ended. Each 
spectrum was fitted in XSpec using a power law, V x, the Galactic 
Njf ( neutral hydrogen column density, taken from iKalberla et al.l 
l2005h and the intrinsic N/f at the redshift provided in Table |2] The 
spectral fits are provided in Table |4] 

The majority of the SGRBs are consistent with having neg- 
ligible intrinsic N// observed in their spectra suggesting they are 
likely to have occured in low density environments. Recently, 
Margu tti et al.! (^012) have compared the distribution of intrinsic 
Ni/ observed in SGRBs to LGRBs finding that they are typically 
consistent with the lower end of the LGRB distribution consistent 
with the higher end of our distribution and we find several candi- 
dates with negligible intrinsic absorption. Some of the sample have 
significant N// values, but it is important to note that detailed ob- 
servations have shown that the optical absorptions found for GRB 
afterglows can be ord ers of magnitude less than that expected fr om 
the X-ray N// values dSchadv et al.ll20ld : ICamDana et al.ll201Ql) . 



3. 3. 2 Optical Afterglows 

A 1 keV observed flux lightcurve showing the prompt. X-ray and 
the most constraining optical observation during the plateau phase 
was created for each burst in the sample. These were produced 
using the simple relation given in equation |7] (assuming a simple 
power law spectrum and a spectral index /3x = r^c — 1) to shift the 
observed fluxes at a measured energy to 1 keV. 
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Figure 10. (a) The 0.3 - 10 keV unabsorbed flux at 100 s versus 1000 s. (b) The 0.3 - 10 keV unabsorbed flux at 100 s versus 10000 s. (c) The 15 - 150 keV 
flunce versus the 0.3 - 10 keV unabsorbed flux at 100 s. (d) The 15 - 150 keV flunce versus the 0.3 - 10 keV unabsorbed flux at 1000 s. Symbols are as in 
Figure |9] 
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Fx was obtained from the time averaged PC mode spectra pro- 
duc ed by the automated anay lsis on the UK Swift Data Centre web- 
site ([Evans^ alJ200l[200^ . The 0.3 - 10 keV observed BAT-XRT 
lightcurves were extrapolated to flux at 1 keV using equation [T] The 
optical magnitudes were converted into flux for the wavelength of 
the optical filter used and then shifted to 1 keV using equation [T] 
As there may be a cooling break inbetwee n the optical and X-ray 
observations (ISari. Piran. & Naravanll 19981) . the two extreme cases 
are taken i.e. /3o = 13 x and f3o = f3x — O.b. The errors on the 
observed optical magnitudes and the errors on r^c are used to de- 
fine the region on the lightcurve that the optical data could reside 
in (dark grey - no cooling break, light grey - cooling break, note 



there is overlap between these two regimes). If the optical and X- 
ray data are consistent, then the X-ray data points should lie within 
the shaded regions for the optical data. 

The 1 keV flux lightcurves for SGRBs fitted with the magnetar 
model are shown in Figure [TT]compared with the most constraining 
optical observation extrapolated to 1 keV. These compare the BAT- 
XRT lightcurve at 1 keV to the most constraining optical obser- 
vation extrapolated to 1 keV. GRBs 051221A, 061201, 080905A, 
080919 and 090426 have optical afterglows which are consistent 
with their X-ray afterglows, but many of these would require the 
most extreme errors on the spectral slope and cooling break. ^55% 
have optical afterglows that are underluminous with respect to their 
X-ray afterglows, signifying either significant optical absorption or 
an extra component in the X-ray afterglow. However, as shown in 
Section 3.3.1 using absorption in the X-ray spectra, the majority of 
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Figure 11. Comparison of the X-ray and optical data for the SGRBs fitted with the magnetar model. These are observed X-ray flux lightcurves at 1 keV with 1 
extrapolated optical observation, light shaded region = optical observation assuming the most extreme cooling break between X-ray and optical and dark shaded 
region = optical observation assuming no cooling break. The references are for the optical observation used. If the X-r ay and optical observa tions are consistent 
with originating from th e same source, the X-ray data points sh ould pass throu gh the shaded regions. GRB 050509B - iBreeveld et akltoOSh - consistent, GRB 
051210 - Jelinek et al." 12005) - inconsistent, GRB 05 1221 A - Soderberg et al. (2006) - optical observations are consistent with X -ray ob servations, GRB 
060313 -fRoming et al. (2006) - inconsistent, GRB 060801 - Brown & Racusin (2006.) - inconsistent and GRB 061201 - .Stratta et"ZI ( l2007h - only consistent 
for with most extreme cooling break and errors. 
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Figure 11 - c ontinued GRB 0707 14A - 1 Chester & Gmpd (l2007h - upper limits inconclusive if there is an extreme cooling brea k, GR B 070724 A - 
Ide Pasquale & Z iaeepour ( 2007) - up per limits inconcl usive if there is an extreme cooling break, GRB 070809 - IChester & Marshall ( l2007h - upper limits 
inconclusive, GRB 080 426 - Gates & ZiaeepouJ (l2008h - inconsistent, GRB 080702 A - .de Pasq uale (2008) - upper limits inconclusive and GRB 080905 A - 
Brown & Paganil ( l2008h - upper limits inconclusive. 
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Figure 11 - continued GRB 080919 -'immler & Ho Uandl l2008h - likely consistent GRB 081024A - Ide Pasquale & Strattal bOOSh - upper limits inconclusi ve 
and GRB 090426 - Gates & Cumminss ( 2009) - optical observations are consistent with X-ray observ ations, GRB 090510 - iKuin & HoverstenI ( l2009l) - 
inconsistent, GRB 090515 - .Siesel & Beardmora (i2009l) - inconsistent and GRB 090621B - ICurranI (l2009h - inconsistent. 
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Figure 11 - continued GRB 09110 9B - lOatesI l2009h - upper limits inconclu sive if there is an extrem e cooling break, GRB 1001 17 A 
- Pasquale, Holland, & OatesI <20ld) - extremely incon sistent, GRB 100625A - [Landsman & Holland (2010h - inconsistent, GRB 10 702A - 
Ide Pasquale & Siesel (2O10h - inconsistent, GRB 101219A - iKuin & GelbordI bOlOh - inconsistent and GRB 110112A - iBreeveld & StamatikosI l201lh - 



mconsistent. 
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Figure 11 - continued GRB 1 1 1020 A - iGuidorzi et alJ bOllh - co nsistent GRB 1 1 1 1 17A - lOates & ManeanoH201lh - upper l imits inconclusive if there is an 
extreme cooling break, GRB 120305 A - [Marshall & Strattal ( I2OI2I) - inconsistent and GRB 120521A - loates & ChesteJ ( 1201 2h - inconsistent. 



the candidates are consistent with occuring in a low density envi- 
ronment. 

In Figure [12] we compare the average X-ray fluxes at 1 keV 
to the optical fluxes extrapolated to 1 keV with (b) and without 
(a) a cooling break in the spectrum. The average X-ray flux was 
calculated during the optical observation. There are several points 
which lie below the black line in both cases, showing there is less 
emission at optical wavelengths than expected. 

To determine if the observed X-ray excess could be caused by 
optical absorption, we compare the optical absorption (Ay) esti- 
mated using the observed X-ray N// to the minimum absorption 
that could explain the difference between the X-ray and optical 
fluxes. The observed spectra during the plateau regime (given in 
Table lU are used when available and the other spectral fits are ob- 
tained from th e automated data products from the UK Swift Science 
Data Centre (Evans et al. We convert the observed X- 

ray N// to optical absorptions using for Milky Way (MW, 

1.8 X 10^^ Predehl&Schmittlll995h. Large Magell anic Cloud 
(LMC, 3.5 X 10^^ : iKoomne ef "l982: Fitzpatrick 1985) and Small 
Magellanic Cloud (SMC, 4.0 x 10^^ ; Martin, Maurice, & Leaueuxl 
Il989l) abundances. Note that there are known to be significant scat- 



ter and uncertainties involve d in this conversion (e.g. ISchadv et alJ 
I2OIQI : ICampana et all 120101) . To obtain the minimum Av^ which 
would be sufficient to explain the difference between the X-ray 
and optical fluxes, the maximum possible optical flux (including er- 
rors and assuming the most extreme cooling break) and the X-ray 
plateau flux are converted to V band magnitudeo The obtained 
optical absorptions are given in Table |5] and plotted in Figure [13] 
Many of the GRBs may be explicable via absorption however we 
note that ^25% of the sample are based on unconstraining optical 
upper limits while some rely on using the most extreme cooling 
breaks and uncertainties. In Figure [13] we also show that if some 
of the host galaxies are more consistent with LMC or SMC abun- 
dances then more of the GRBs cannot be explained via absorption. 
Results obtained by Schady et al. (2010) for LGRBs, also suggest 
that may be an order of magnitude higher for GRB host galax- 
ies, in which case even more GRBs in the sample would not be ex- 
plicable via absorption. Despite all the uncertainties involved in this 
calculation, 8 GRBs in the sample clearly cannot have the differ- 
ence between their X-ray and optical fluxes explained via absorp- 

^ using the webtool: [littp://www.stsci.edu/hst/nicmos/tools/conversion_form.html| 
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Figure 12. The optical flux shifted to 1 keV is plotted against the average X-ray observed flux during the optical observation also shifted to 1 keV. The solid 
black line represents where these are equal, as expected if they are consistent with each other. In (a) we assume there is no cooling break between the optical 
and X-ray observations and in (b) we assume the most extreme cooling break. Symbols are as in Figure |9] 




Ay predicted using predicted using A^ predicted using 

Figure 13. A plot comparing the minimum optical absorption, required to explain the difference between the X-ray and optical absorptions to those 
predicted using the X-ray ^h- Unless the optical data are already consistent with the X-ray observations, all data points are lower limits given the assumptions 
made. These plots are for three different abundances: (a) Milky Way, (b) Large Magellanic Cloud and (c) Small Magellanic Cloud. Data points lying above 
the black line cannot be explained by simply using optical absorption. Symbols are as in Figure |9] 



tion(GRBs 060313, 061201, 090510, 090515, 1001 17A, 100625A, 
100702Aandll0112A). 



This analysis shows that at least some of the GRBs in this 
sample are consistent with there having an additional X-ray com- 
ponent. This may provide supporting evidence of energy injection, 
although energy injecti on is thought to cause an increase in flux 
at all wavelengths (e.g. ISari & Meszarosll2000l however this also 
depends upon the electron energy distribution). 



Although there is some evidence that the magnetar candidates 
have additional X-ray emission, it is not known what spectrum is 
expected from a newly formed magnetar and hence we cannot com- 
pletely discount those whose optical emission is consistent with 
their X-ray emission. 



4 DISCUSSION 

4.1 The sample of SGRBs 

Here we discuss some particular SGRBs and then the sample as a 
whole. 

GRB 070809 is one of the best fitting stable magnetar can- 
didates and lies within the allowed regions. This GRB had a faint 
optical afterglow an d is offset by 20 kpc from a galaxy at z = 0.219 
dPerlev et alJl2008h . making it an ideal candidate for a magnetar 
formed via the merger of two NSs. However it is important to be 
cautious about this candidate host galaxy association as the likeli- 
hood that this is an unrelated field galaxy is 5 - 10% (Tunnicliffe et 
al. in prep). 

GRB 061201, with a spin period of ^^16 ms, fits the magne- 
tar model well but is spinning slower than expected. However the 
redshift used relies on the correct host galax y identification which 
remains highly uncertain (IStratta et aLll2007l ; Tunnicliffe et al. in 
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Unstable 


081024 


Yes 


7 


Yes 


Unstable 
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Unstable 


090621B 


7 


Yes 
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No 


Stable 


1001 17A 


Yes 


Yes 


No 
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100625A 
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Yes 
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100702A 


Yes 


Yes 


Yes 


Unstable 


101219A 


Yes 


Yes 


Yes 


Unstable 


110112A 


7 


Yes 
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Stable 


111020A 


Yes 


No 


No 


Stable 


111117A 
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Stable 


120305A 
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Yes 


Unstable 


120521A 


Yes 


Yes 


Yes 
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Table 6. A summary showing the main features studied. This gives best magnetar candidates found and the possible candidates. "Expected region" : fits within 
the required parameter space in Figure |9](? = could fit with various assumptions), "Extra component" : there is evidence of an extra component in the X-ray 
afterglow which is not observed in the optical note this could also be due to absorption (? = borderline case or optical upper limit not constraining), "Pred icted 
region " : do the values for the plateau luminosity and the plateau duration, calculated using equations [T] and |2] lie within the predicted region in Metzge r et al.l 
i201lh ? (? = outside region but would fit with reasonable assumptions). "Stable/Unstable" : whether the magnetar is stable or if it collapses to form a BH (? = 
would be fitted well by either case.) 



prep). If it actually occured at a higher redshift than used in this 
analysis it would lie within the expect ed region. Additionally, the 
approximate 10 ms limit imposed bv Usovl(ll992) is dependent on 
the initial radius of the collapsing object and the radius of the fi- 
nal NS. This limit is also derived for the model involving AIC of 
a WD. Therefore there is some level of flexibility in this imposed 
limit. We still consider this, and other GRBs close to this boundary, 
to be potential candidate magnetars. 

GRB 05 1221 A is consist ent with having energy injection in its 
lightc u rve out to '- ^ 2 x 1 0^ s (iBurrows et al l2006l : ISod erberg et al. 
l2006h . iFan &3a] (l2006h explained this as energ y injectio n from 
a magnetar. Our model fits this GRB very well. iJin et al 
proposed an alternative two jet model to explain the lightcurves 
without requiring additional energy injection. 

GRB 060313 has been included in the magnetar sample by 
ignoring the first 50 - 200 s of the lightcurve due to the flar- 
ing activity, this gives a good fit to the later data but this result 
should be treated with caution. Flares could be associated with 
on-going accretion on to the newly formed magnetar. Alternatively, 
iDai et alJ (l2006h and iGao & FanI (l2006h suggest that the X-ray 
flares originate from reconnection of twisted magnetic fields within 
the NS. Margutti et al. (2011) have conducted a systematic study 
into SGRB flares, including the flares observed in GRB 060313, 
and concluded that the flares are consistent with a central engine 
origin. 

Included in this sample are SGRBs whose progenitors are 



subject to significant debate, particularly GRB 090426 at z^2.6 
which c ould have orginated f r om a coUapsar instead of a binary 
merger (lAntonelli et all l2009l: iLevesque et~aD l2010l : iThone et al] 
l201ll : IXin et aLll201lh . GRB 090426 fits the model well, irrespec- 
tive of the progenitor, but the progenitor debate is important to note 
as we are specifically studying possible NS binary merger progen- 
itors. 



Interestingly, 12 out of the 28 magnetar candidates require col- 
lapse to a BH. This implies that, if these SGRBs are making magne- 
tars, they only collapse to a BH in a small number of cases. Com- 
paring the derived plateau durations and the collapse times pro- 
vided in Table [3l the magnetar typically (but not always) collapses 
to a BH after the plateau phase, i.e. when the magnetar has spun 
down significantly. The only exception to this is GRB 101219A 
where collapse occurs prior to the end of the plateau phase, how- 
ever the collapse time and end of the plateau are consistent within 
errors. The collapse time is related to the mass of the magnetar and 
the spin period at which the differential rotation can no longer sup- 
port gravitational collapse. The discrepancy between collapse time 
and plateau duration are hence likely to be reliant upon the mass of 
the magnetar. Additionally, there may be ongoing accretion on to 
the magnetar (remnants of the merger) which may raise the mass of 
the magnetar above the critical point prior to significant spin down. 
Interestingly, those candidates which collapse to form a BH and 
are within the allowed (unshaded) region of Figure |9]have a higher 
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GRB Vx Galactic N// Restframe Intrinsic N// 

(1020 cm-2) (1020 cm-2) 



Magnetar candidates 



051221A 


2Q4+0.14 


5.70±0.37 




060313 


1 ^1 -1-0 16 
^•"^-0.13 


5.00±1.17 


'^•'^'^-0.00 


060801 




1.40±0.31 


OQ Q+68.8 
^^•^-29.9 


070809 


1 70+0. 83 
^•'■^-0.43 


6.40±0.17 


9 QC + 14.9 
^•^-^-2.95 


080426 


. Qo+0.29 
^•^-^-0.27 


37.0±4.19 




080919 


9 90 + 1.02 
^•^■^-0.84 


26.0±3.78 




090426 


^•^■^-0.11 


1.50±0.11 


00+^6 
'^•'^'^-0.00 


090510 


^•-^^-0.19 


1.70±0.11 


10 0+160 


090515 


^•^^-0.24 


1.90±0.25 


13 1 + ^^-^ 
^■^•^-10.5 


101219A 


1 65+°-^2 


4.90±0.87 


+ 26.7 

^0.8_2Q 4 


111020A 


2 5^+1.69 


6.89±0.48 


7 94+7.90 
'•^^-7.90 


120305A 


1 94+0.21 
^•^^-0.20 


11.3±0.70 




120521A 


1 61+°-^^ 


20.80±1.69 


1 9 + 14.2 
^•^-1.2 



Possible candidates 



050509B 


1 92+1-09 

^•^^-0.60 


1.60±0.04 


8-OOt«3:J°o 


061201 


1 44+0.20 

^•^^-0.19 


5.20±1.58 




0707 14A 


.^+0.37 
^•^^-0.35 


9.20±1.25 


2141^1;? 


080702A 


. ^-7+0. 85 
^••^'-0.76 


15.0±1.50 


19^- + 251 
i^J_121 


090621B 


9 3Q+1.60 


19.0±1.96 


42 8+^°^ 


091109B 


1 96+0-64 
^•^"-0.43 


9.20±0.96 


.43 + 27.9 
^^•-^-14.5 


110112A 


2 07+°-^^ 


5.50±0.40 


7 oz,+ 12.7 


111117A 


9 .0+0.39 
^•^■^-0.36 


3.70±0.15 


39.8l«?:^ 



Table 4. The 0.3 - 10 keV spectral fits for the derived plateau durations 
given in Table [3] These are the SGRBs in the magnetar sample which 
have X-ray data during the plateau phase. Provided are the photon index, 
^ X, plateaus the Galactic and the restframe intrinsic using the 
redshifts provided in Table |2] 



magnetic field for a given spin period than the candidates which do 
not collapse to a BH. 

Many of the magnetar candidates lie within, or near to, the 
predicted plateau lu minosity and duration regions for newly formed 
magnetars given in lMetzger et alj (|201 ih when considering uncer- 
tainties due to redshift, efficiency and beaming. However, there are 
candidates whose plateaus are significantly shorter than predicted 
or at a lower luminosity. Our analysis and that of iMetzger et al.l 
(l201lh assumes a NS mass of 1.4M0 and this is likely to be signif- 
icantly higher for a NS merger progenitor (e.g. 2.IM0). This has 
a small affect on the values of the magnetic field strength and the 
spin p eriod calculated in our model (as shown in Rowlin son et al.l 
l2010ah but does not significantly af fect the predicted reg ions for 
plateau luminosity and duration from^ Metzger et alj (|201 lb . 

A summary of the properties of the whole magnetar sample 
are shown in Table [6] 



4.2 Accretion Effects 

In our analysis we have not accounted for any ongoing accretion 
onto the magnetar from the surrounding torus of material formed 
during the merger. This could significantly effect the results ob- 
tained, especially if accretion increases the NS mass to more than 
can be supported as this results in collapse to a BH. Addition- 
ally, accretion could explain flares observed overlaying the plateau 



GRB Minimum Ay MWAy LMC Ay SMC Ay 



Magnetar candidates 


051221A 


0.00 


1 32+0-42 

— 0.39 


0.68+0-21 

— 0.20 


0.59+0-i? 

— 0.17 


060313 


2.15 


0.28+°-39 

— 0.07 


0.14+°-20 

— 0.03 


0.13+°18 

— 0.03 


060801 


0.91 


1.75+3-86 

— 0.85 


0.89+?,-9^ 

— 0.43 


0.78+000 

— 0.38 


070724 


0.32 


0.74+?,-^2 

— 0.73 


0.38+0-^^ 

— 0.37 


0.33+0-^0 


070809 


0.00 


0.52+°-?4 

— 0.17 


0.27+°-*3 

— 0.09 


0.23+°-i5« 

— 0.08 


080426 


1.10 


3.85+2-00 

— 1.66 


1.97+1-02 

— 0.85 


1.73+0-89 

— 0.74 


080905A 


0.46 


1 61+0-95 


0.83lg;« 


72+0-43 

^- '^-0.33 


080919 


1.12 


7 -29 + 7.25 
/"^^_4.45 


3 74+3.71 

— 2 27 


00+3. 24 
-'•^O-l 99 


081024 


0.00 


7 09+4.47 
'•-^^-3.02 


3 74+2:29 
■^•'^-1.54 


3.28+2-00 


090426 


0.00 


o.o8i§:°j 


QQ4+1.03 

'^•'^^-o.oo 


QQ4+o:9o 

'^•'^^-o.oo 


090510 


2.12 




Q 33+0.46 
'^•■^-^-0.29 


0-29tS:f5 


090515 


3.55 




Q 43+0.34 

'^•^■^-0.31 


0-38l°o:f7 


1001 17A 


6.98 


9 00 + 1.96 
^•■^■^-1.45 


1 IQ+I.OO 
^•^^-0.74 


1 04+0-88 
^•^^-0.65 


100702A 


7.17 


4 13 + 2.07 
^•^■^-1.79 


9 1 1+1.06 
^•^^-0.91 


1 o^r + 0.93 
^•"^-0.80 


101219A 


2.30 


3.4511:- 


^•''^-0.61 


1 54+0-69 

^••^^-0.53 


111020A 


0.00 




Q 49+0.24 
^•^^-0.24 


Q 37+0.21 
'^•■^'-0.21 


120305A 


3.47 


^^•'^-1.49 


3 44+0.93 
■^•^^-0.76 


3 01+°-82 
.^•'-'^_0.67 


120521A 


0.78 


1 90+0.89 

^•^■^-0.16 


Q 63+0.45 
'^•'^■^-0.08 


55+°-^° 

'^•'^^-0.07 


Possible candidates 


050509B 


0.00 


Q 34+0.45 
'^•-'^-0.45 


n 97+0.23 
'^•^'-0.23 


Q 94 + 0.20 

'^•^^-0.20 


051210 


0.00 


2 45+1.01 

^•^-^-1.01 


1 9^-+0.51 
^•^-^-0.51 


1 10+0-45 
^•^^-0.45 


061201 


2.19 


Q ^7+0.33 
'^•^ -0.31 


Q 34+0.17 
^•■^^-0.16 


Q3Q+0.15 

'^•■^'^-0.14 


0707 14A 


0.00 


12 47+2-96 


6.3811-J2 


^- ^-0 + 1.33 
j.Jo_-|^ 27 


080702A 


0.00 


7 09 + 14.11 


4 00+'^-2i 

^•'^'^-3.50 


3 50+6-31 

■^••^'^-3.06 


090621B 


0.99 


3 43+6.14 

■^•^^-2.50 


1 77+3.14 
^•^ ^-1.28 


1 55+2-'^5 


091109B 


0.00 


1 09 + 1. 61 
^•■^^-0.86 


n AQ + 0.82 

0.68^044 


Q3Q+0.72 
^^•^^-0.39 


100625A 


3.49 


12+0-42 

'^•^^-o.oo 


0.06ig:gJ 


05+0-19 


110112A 


1.29 


Q 73+0.73 
'^•''^-0.46 


0.381-.37 


(\ 00 + 0.33 
^•■^■^-0.21 


111117A 


0.00 


940+3.90 
^•^■^-1.76 


1 94+2.00 

^•^^-0.90 


1 09+1-'^^ 
^•^^-0.79 



Table 5. The minimum optical absorption, Ay , is the absorption required 
for the optical observations to just be consistent with the X-ray observations 
(0 means they are already consistent). The MW (Milky Way), LMC (Large 
Magellanic Cloud) and SMC (Small Magellanic Cloud) absorptions are the 
values predicted using the X-ray • 



model. Flares may also be associat ed with ongoing magnetar activ- 
ity as described in Da i et aP (l2006h . 

Piro & Ott (201ll) studied the affect of accretion onto mag- 
netars formed during SNe, however their results are also applica- 
ble to magnetars produced from neutron star binary mergers. The 
main difference for mergers is the significantly reduced reservoir 
of material available for accretion and the different accretion rate. 
In this section , we as sume the simplest accretion rate published by 
iMetzger et al.l (|201Q|) assuming that accretion starts at 0.16 s af- 
ter the trigger time, giving a total accretion disk mass of ^ 0.3 
M0 . Accretion onto the magnetar occu rs when the the pr opeller 
regime ends, given by equation [8] from IPiro & OttI fcoill) where 

M < 6.0 X 1{)-''^iI^M~1^''p~]/^Mqs-^ (8) 

As before, we assume an initial NS mass of 1.4 M© and radius 
of 10^ cm. In Figure [T4k we show the accretion rate as a function 
of time after formation. In Figure [T4l ) we show the evolution of 
the spin period of two different magnetars (using the parameters 
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100 
Time (s) 

Figure 14. (a) The accretion rate as a functi on of time ass uming the accre- 
tion rate for a compact binary merger Metzg er et al 1 (l2010h starting at 0.16 s 
after the trigger time giving a total accretion disk mass of ~0.3 Mq . (b) The 
evolution of the spin period of the magnetar for the two accretion rates, red 
- the magnetar predicted for GRB 060313 and blue - GRB 090515. Solid 
lines include accretion and dashed lines have no accretion. In these plots, 
accretion has a very small or negligible effect, (c) The amount of rotational 
energy available in the magnetar for each case. 



for GRBs 060313 and 090515 as these have contrasting magne- 
tar properties) assuming there is either accretion onto the magne- 
tar or no accretion. When there is significant accretion (e.g. GRB 
090515) it can marginally prevent spin down and increase the rota- 
tional energy (Figure [T4b ) available, although these are negligible 
effects for the low accretion rates considered. 

It is worth noting that accretion would potentially have a very 
large effect on the results obtained fo r LGRB magnetar candidates 
(e.g. the sample in Lyons et al.ll201Ql) as these are thought to have a 
significantly higher mass accretion disk and an accretion rate sim- 
ilar to that proposed by Piro & Ott (2011). In that case, the energy 
reservoir could reach values in excess of 10^^ ergs for particu- 
lar combinations of the initial conditions. This additional energy 
source could be a potential explanation for large flares observed in 
some of the LGRB candidate lightcurves (e.g. Margutti et al. 201 1). 

In Figure [15] we show the total mass accreted after the pro- 
peller regime has ended. The linear correlation between the dura- 
tion of the propeller regime and the mass accreted is caused by the 
relationship: M oc t~^^^ (i.e. the sooner the propeller regime ends, 
the greater the mass that can be accreted). The candidates which 
accrete the most mass are those which also collapse to form a BH 
within a few hundred seconds, leading to the suggestion that accre- 
tion is an alternative to drive this collapse. Typically, the magnetar 
is thought to collapse when the fast rotation can no longer support 
the mass of the magnetar. The stable magnetar outliers are GRBs 
100625 A and 1001 17A which were also well fit by the unstable 
magnetar model but we chose the stable model to reduce the num- 
ber of free parameters. Additionally GRB 090426 is again a clearly 
stable magnetar candidate which is separate from the other stable 
candidates. 
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Figure 15. The amount of mass accreted by the magnetar against the dura- 
tion of the propeller regime. The dashed line represents the maximum mass 
available in the accretion disk and is 0.3 an upper limit for the amount 
of mass which can be accreted. Symbols are as in Figure |9] 




10 100 
E,o(10'' erg) 

Figure 16. The energy emitted during the plateau phase, calculated using 
the fits in Table [3] compared to the isotropic energy emitted during the 
prompt phase (1-10000 keV). Symbols are as in Figure |9] 



4.3 Energy Constraints 

Including all of the possible candidates, the SGRBs in our sam- 
ple can be fitted with the magnetar model. In Table [3] we show 
the isotropic energy released during the prompt emission phase of 
the GRB. These values tend to be consistent with the maximum 
expected energy outpu t from the magneta r central engine model, 
Eiso < 3 X 10^^ erg (iMetzger et allboili) . Within the uncertain- 
ties many of the magnetar candidates are consistent with this limit 
while some others exceed it. However, we have not corrected for 
beaming and had to assume redshifts in many cases. Not correct- 
ing for beaming will undoubtedly affect these results by increas- 
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ing the spin p eriod and the magnetic field strengths as shown in 
iRowlinson et a l. ( 2010a). Beaming, with a half-opening angle of 
30° , has been shown to form via the formation of an ordered ma- 
gentic field during the merger of two 1.5 M© NSs which collapse to 
form a BH (iRezzolla et al.ll201lb . However, the beaming angles of 
SGRBs and associated magnetars remain unconstrained (see recent 
work on SGRB jets by Fong et al. 2012). With a reasonable beam- 
ing correction, all of the GRBs which exceed the energy constraint 
would lie well below the maximum expected energy output. 

Another consideration is that Eiso oc Mi.4Po7-3' mag- 
netars can have masses up to 2.IM0 then the maximum energy 
output could be as high as E^so ~ 1 x 10^^ erg. 

In Figure [161 we show the energy emitted during the magnetar 
plateau phase (the plateau luminosity multiplied by the duration 
from Table [3] these values were calculated from the fitted B15 and 
P_3 using Equations |2] and [B against the isotropic energy emitted 
during the prompt emission. Only five GRBs which fit the magnetar 
model emit more energy during the plateau phase, GRBs 051210, 
070724A, 070809, 090515 and 100702A. 

We have also assumed 100% efficiency in the conversion of 
rotational energy into EM radiation. This will not be the case and 
assuming a lower efficiency would act counter to the effect of 
any beaming, in the sense of reducing the inferred spin period 
and the magnetic field strengths. For example, GRB 090515 has 
B - 1.4 X 10^^ G and P - 2.3 ms assuming 100% efficiency, 
at 10% efficiency these drop to P - 4.4 x 10^^ G and P - 0.73 
ms. Given the uncertainties in both beaming and efficiency, we note 
that the real values of the magnetic field strength and the spin pe- 
riod may be uncertain by at least a factor of 3. 



4.4 Gravitational Wave Signals 

Systems of the kind we have considered represent interesting 
sources of gravitational waves as there are predicted signals for all 
of the stages this system can go through: inspiral, magnetar spin- 
down and final collapse to BH. In Table [T] we show the distances 
out to which each phase would be visible, assuming the ampli- 
tude Qi) of the gravitational waves is inversely proportional to dis- 
tance for Advanced LIGO (AdLIGO, with a sensitivity of /i 4 x 
10~^^), the Large Cryogenic Gravi tational Telescope (L CGT, com- 
parab le sensitivity to AdLIGO; iKuroda & LCGT C ollaboration 
20101) and the Einstein Telescope (ET, /i - 3 x 10"^^ Hild et al. 
201 ll) . The gravitational wave amplitude is quoted for a distance 
of z ~ 0.1 or 390 Mpc. The magnetar phase prediction is an up- 
per limit assuming a spin period of 1 ms, /45 = 1-5 for a binary 
merger progenitor, and an ellipticity e — 1. AdLIGO predictions 
byLAbadie et aL| (I201QI) are for NS-NS mergers. 

Using the lowest and maximum possible rate s for NS-NS 
merge rs per Milky Way Equivalent Galaxy from lAbadie et al.l 
feOlOh . it is possible to predict the number of unstable magne- 
tars (i.e. one source giving 2 distinct gravitational wave signals) 
we might expect to detect with AdLIGO and ET. To detect all 
the stages for the formation and collapse of a magnetar, AdLIGO 
would require it to be at a distance ^100 Mpc and ET would re- 
quire ^^1300 Mpc. Within these volumes there is predicted to be 
a NS-NS merger rate of 2 x 10"^ - 0.08 yr"^ for AdLIGO and 
10 - 4 X 10^ yr~^ for ET. However, the rates need modification as 
not all NS-NS mergers will lead to an unstable magnetar which will 
give both signals. From the analysis in this paper, only 1 1 SGRBs in 
the total sample of 28 magnetar candidates (39%, assuming NS-NS 
mergers always produce a magnetar) are thought to form unstable 
magnetars, giving rates of 8 x 10~^ - 0.03 yr~^ for AdLIGO and 



4 - 2 X 10^ yr"^ for ET. Therefore, it is unlikely that AdLIGO or 
LCGT will observe both the formation and collapse of an unstable 
magnetar b ut ET should de tect many cases. On a more optimistic 
note, Bauswein et al.l (l2012h estimate that AdLIGO will be able to 
detect a post-merger signal associated with a newly formed massive 
NS with a rate of 0.015 - 1 2 yx~ \ 

IShibata & Taniguchil (l2006h also study different masses rel- 
ative to the maximum mass of a NS. They determined that if 
M < Mmax then the NS will emit gravitational waves during the 
magnetar spindown phase until it is a stable sphere and collapse 
to a black hole is dependant on the gravitational wave emission 
(possibly collapsing within 50 ms) or on forces such as magnetic 
breaking. In this case, they predict that advanced gravitational wave 
detectors will be able to observe these events out to 50 Mpc using 
detectors such as AdLIGO. Alternatively if M ^ Mmax, then it 
collapses rapidly to spherical shape and hence is more likely to 
create a stable NS which may collapse at late times due to mag- 
netic breaking. The gravitational waves from the more massive NS 
would be detecta ble to 10 Mpc. 

In both [Baiotti. Giacomazzo. & RezzoUal <2008l) and 
IShibata & Taniguchil ( l2006h . instabilities in the NS formed 
by a compact merger produce detectable gra vitational waves i n 
contrast to the sphe r ical c ollapse model of IPiro & Ottl (1201 ih . 
However IPiro & Ot3 (l201lh showed that accretion may have an 
important affect on the gravitational wave signal. Therefore, these 
objects are potentially important sources of gravitational waves 
and further analysis combining all these factors and the new limits 
on maximum NS masses is required. 

The predictions by Metzg er et al ](|20T 1) do not take into ac- 
count the loss of energy via gravitational waves and this may play 
a significant role for the formation of a magnetar via the merger 
of two NSs. Som e of our candid ates have shorter plateau durations 
than predicted bv lMetzger et al. l ( 201 1 ) however if the energy losses 
via gravitational waves are more significant then the magnetar will 
spin down more rapidly. 



5 CONCLUSIONS 

We have analysed the BAT-XRT lightcurves of all the Swift GRBs 
with prompt durations Tqq ^ 2 s detected until May 2012. About 
half of these SGRBs require fitting with a broken powerlaw model 
showing a plateau phase. Although the plateau phases show many 
similarities with those observed in LGRB lightcurves, they are typ- 
ically orders of magnitude earlier. The initial temporal indicies (ai 
and a2) are comparable to those found for the "canonical" LGRBs 
but there is much more variation in the final decay (as). The cor- 
relation between luminosity and duration of the plateau phase is 
found to be consistent with the identified correlation for "canoni- 
cal" LGRB lightcurves identified by Dainotti et al. ( 2010). 

Following on from the study of GRB 090515, this work has 
shown that the X-ray lightcurves of some SGRBs considered could 
be explained with energy injection from a magnetar which can col- 
lapse to form a BH. 18 firm candidates (64%) and 10 possible can- 
didates were found. Of the 18 firm candidates, 10 are thought to 
collapse to form a BH and when including possible candidates, 1 1 
out of 28 magnetar candidates may collapse to form a BH. This 
implies that 29-56% of events forming magnetars would collapse 
to a BH within the first few hundred seconds. In some cases the 
magnetar plateau phase is not directly observed as it occurs prior 
to the XRT observations. This predicts plateau emission that may 
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Phase Citation Predicted Amplitude Distance used AdLIGO/LCGT limit ET limit Amplitude at z~0. 1 

(h) (Mpc) (Mpc) (Mpc) (h) 

Inspiral Abadie et al. (2010) 4 x IQ-^^ 445 445 5900 4.6 x 10-^4 

Magnetar Spindown Com & Meszaros (2009) 7 x 10-^4 iqo 175 2300 1.8 x 10-^4 

Collapse to BH Novak (1998) 4 x IQ-^s iq loO 1300 1 x 10-^4 

Table 7. Gravitational wave predictions for the three different regimes in this magnetar model and applied to futu re observatories. The distances quoted are 

luminosity distances. The magnetar spindown values are calculated using Equation 14 in lCorsi & Meszaro'sl ( l2009h . 



be observable with future missions that are able to slew faster than 
Swift. 

The X-ray fluxes at 1000 s and 10000 s are typically higher 
for the stable magnetar candidates. The late time fluxes are signifi- 
cantly lower for the unstable magnetar cases. There is excess emis- 
sion in the X-ray afterglows not observed in the optical afterglows 
for many of the magnetar sample. Many of the magnetar candidates 
lie within or close to the predicted regions for plat eau luminosity 
and du ration for newly formed magnetars given in iMetzger et alj 
(I2OIII) . 

Accretion onto the newly formed magnetar formed by a NS- 
NS binary merger has a negligible affect on the spin periods and 
hence the rotational energy budget of the magnetar. However, it 
can be shown that accretion can have a significant affect for collap- 
sar progenitors. This may explain late time flares for collapsar pro- 
genitors and our calculations suggest the rotational energy budget 
could exceed 10^^ erg for some combinations of initial spin peri- 
ods and magnetic fields. The unstable magnetar candidates, those 
which collapse to form a BH, are potentially accreting more mate- 
rial than the stable candidates. We suggest this is an additional solu- 
tion for why they collapse at late times which would work alongside 
the theory that the magnetar spins down to a critical point where it 
can no longer support its mass using rotation. 

These objects are highly interesting targets for future gravita- 
tional wave observatories as they are predicted to emit gravitational 
waves during merger, the magnetar phase (likely to be increased via 
accretion and bar mode instabilities) and, in some cases, the final 
collapse to form a BH. In this paper, we have focused on NS-NS 
merger progenitors, however the accretion induced collapse of a 
WD could also produce a SGRB and leave behind a rapidly rotat- 
ing magnetar with similar X-ray emission properties. Among other 
observational signatures, the very different gravitational wave sig- 
nals between these events may someday allow these progenitors to 
be distinguished however the inspiral remains the most luminous 
phase of gravitational wave emission. 

For the candi d ates which form a stable magnetar, 
iDuncan & ThompsonI ^1992) showed that the amount of en- 
ergy available for an SGR giant flare E oc 3 x 10^^ B^^ erg. 
Hence a young magnetar, with a magnetic field of B15 ^ 10, 
could produce a giant flare with an energy of 3 x IC^^ erg. This 
value is comparable to the isotropic energy of some SGRBs (e.g. 
GRB 080905A at ^ - 0.12, R owlinson et al.l ii2010b) so would 
be observable in the local universe. Both of the merger and giant 
flare events are very rare, however considering these models it is 
possible (although very unlikely) that in the future we may have 
two spatially co-incident SGRBs. This has also been proposed for 
LGRBs b\ lGianniosl (12010) and they suggest that these magnetar 
candidates could be identified by discovering an old spatially 
coincident radio GRB afterglow in nearby galaxies. 

We have shown that a model of SGRB production from binary 
NS mergers that result in the formation of a magnetar can explain 
the plateaus seen in many SGRB X-ray lightcurves. Although this 



is not conclusive proof of such a model, it would tie in to the evi- 
dence for late time central engine activity in SGRBs and may have 
important observational consequences. 
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